Geometry-induced reduction of the critical current in superconducting nanowires 
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Reduction of the critical current in narrow superconducting NbN lines with sharp and rounded 
bends with respect to the critical current in straight lines was studied at difTerent temperatures. We 
compare our experimental results with the reduction expected in the framework of the London model 
and the Ginsburg-Landau model. We have experimentally found that the reduction is significantly 
less than either model predicts. We also show that in our NbN lines the bends mostly contribute 
to the reduction of the critical current at temperatures well below the superconducting transition 
temperature. 



I. INTRODUCTION 



The experimental realization of the depairing critical 
current in small superconducting specimens was a sub- 
ject of intense efforts since the beginning of the 60's. The 
problem was solved by showing that in a one dimensional 
wire or strip the depairing critical current can be reached 
when certain experimental care is takso^—. Specifically, 
in a thin strip with a width less than the magnetic pen- 
etration depth, the depairing current may be reached if 
magnetic vortices cannot enter the strip. Although a 
potential barrier exists that hampers permanent vortex 
entrjii, there is a non-zero thermodynamic probability 
that a vortex will be thermally excited over the barrier 
and appears in the strip'*. In this case the Lorentz force 
caused by the current will steer the vortex to either side 
of the strip. Strong enough pinning may freeze vortices 
in the strip until the current becomes sufficiently large in 
order to tear them off from the pinning centers. In case 
of sufficiently strong pinning the critical state can still 
be caused by depairing. However the measured value of 
the critical current will be lower than the theory predicts 
due to a decrease of the superconducting cross-section 
by normal cores of the pinned vortices. In narrow strips 
with widths less than 4.4^, where ^ is the superconduct- 
ing coherence length, vortices cannot exist^. In this case 
the depairing critical current can be realized regardless 
of the pinning centers in the specimen. The magnitude of 
the depairing current depends on the purity of the super- 
conductor and the degree of localization. For extremely 
pure or dirty superconductors in the local limit analyti- 
cal expressions for the critical current near the supercon- 
ducting transition can be derived, while for intermediate 
cases as well as for the whole temperature range only 
numerical evaluations provide the correct values^. 

For a variety of superconducting electronic applica- 
tions, the highest achievable current is an important as- 
pect. For example, the detection efficiency of a super- 
conducting nanowire single photon detector (SNSPD)- 
grows when the operation current increasesi^. There- 



fore, knowledge of the nature of the experimentally mea- 
sured critical current is essential for device optimisa- 
tion. Although the two-fluid type of the temperature 
dependence of the critical current was often observed 
in earlier experimental^, this fact alone does not jus- 
tify that the depairing current has been realized. Even 
if the line width remains unchanged, any deviation from 
the straight-line geometry further reduces the measurable 
critical current due to current crowding and correspond- 
ing local increase in the current density near bends or 
curves. In order to fill an area larger than the optical 
wavelength, the nanowire in SNSPD is usually patterned 
as a meander with sharp 180° turns. The thickness d and 
the width w of typical SNSPD nanowires satisfy the con- 
ditions d < ^ and w ^ Acff where AcH' — X^d^^ and A is 
the magnetic penetration depth. Therefore the current 
density remains uniform over the wire cross-section in 
the straight portions of the meandering nanowire. The 
uniformity is disrupted at the turns where the current 
crowds at the inner edges. The crowding increases the 
local current density above the mean current density in 
the straight portions of the wire. Therefore, the poten- 
tial barrier for the vortex entry first disappears near the 
turns and this decreases the measurable critical current 
in comparison to the straight wire of the same cross- 
section. Another consequence of the current crowding 
near turns is that the practically achievable ratio of the 
operation current to the depairing current in the straight 
portions remains less than it could be in a wire with- 
out turns. This makes straight portions less effective in 
detecting photons at wavelengths where vortex-assisted 
detection^^ takes place. Hence, decreasing the strength 
of current crowding near turns would greatly improve the 
spectral range of these detectors. 

The influence of turns with various shapes on the crit- 
ical current has been studied theoretically in the frame- 
work of the London equations for the transport current 
and the screening current around a vortex^^. Another 
approach, which we also involve for comparison here, 
is the numerical solution of the Ginsburg-Landau (GL) 
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equatioiiis for the superconducting line with turns. This 
allows one to find the current at which the superconduct- 
ing state becomes unstable. We compare predictions of 
these two approaches with the experimentally observed 
reduction of the critical current in superconducting lines 
with different bends. We show that in qualitative agree- 
ment with both theories the reduction increases with the 
increase in the bending angle and the decrease in the 
arc radius of the inner edge of the bend. However, the 
amount of the reduction and its temperature dependence 
noticeably differ from theoretical predictions. 



II. EXPERIMENT 

We have studied the effect of bends on the critical 
current of superconducting NbN lines which had nom- 
inal width of 300 nm and a thickness close to 10 nm. In 
SNSPDs from NbN films, the widths and thicknesses of 
nanowires are typically less than 100 and 5 nm, respec- 
tively. The enhanced cross-section of our lines greatly 
improved reproducibility of fabrication while it still satis- 
fied typical SNSPD conditions d < ^ and w ^ Aoff. Also, 
the higher critical currents significantly simplified mea- 
surements and improved experimental accuracy. More- 
over, since the amount of the reduction in the critical 
current due to bends is lager for structures with larger 
w/^ ratio, an increased width allows for easier observa- 
tion and analysis of this effect. The samples had layouts 
as shown in insets in Figs. [T] and [H They included two 
symmetric bends connected by portions of straight lines 
to each other and to enlarged contact pads. Tapered 
connections to contact pads decreased current crowding 
near the contacts. Samples from the batch #1 consisted 
of lines with two sharp bends whose bending angle 9 var- 
ied from 9 = 90° to 6* = 0°; the latter angle corresponds 
to a straight line. Batch #2 consisted of lines with two 
rounded 90° bends. Rounding was made by arc segments 
with different inner radii r, which varied from just a few 
to 500 nm, and the outer radii r + w, correspondingly. In 
order to reduce the dispersion of the experimental data, 
we prepared and tested several (typically 3 to 4) nomi- 
nally identical samples with each layout. To analyse the 
effect of bending angle and radius at 4.2 K we used the 
critical current densities averaged over all samples with 
the identical layout. 

The resistive superconducting transition was measured 
for each sample at a bias current less than 1 pA. We 
defined the transition temperature Tc as the tempera- 
ture corresponding to the 0.1% of the normal state resis- 
tance at 20 K. The critical currents were measured in the 
temperature range from Tc to 4.2 K in the current-bias 
mode. We used slow, a few second long current sweeps 
and battery driven electronics with reduced noise level. 
This regime mimics the operation of SNSPDs and is rel- 
evant if results are to be used for further detector devel- 
opment. However, it is subject to electromagnetic inter- 
ference and fast fluctuations whose effect reduces the ac- 



curacy of determining the critical current near the tran- 
sition temperature. At temperatures T > 0.95Tc, where 
the current-voltage characteristics were non-hysteretic, 
the critical current Ic was determined using a 100 jiV 
criterium emerging from the voltage resolution of our ex- 
perimental setup. At lower temperatures Ic was defined 
as the bias current at which bridges showed clear jump 
from the superconducting to the normal state. 



A. Technology 

The films were deposited on R-plane optically polished 
sapphire substrates by reactive magnetron sputtering of 
a pure Nb target in a gas mixture of Argon and Nitro- 
gen. The substrates were kept at 750° C during film de- 
position. Details on the NbN thin film deposition and 
characterization of the films can be found in ref. 1151 . The 
lines were defined by electron beam lithography and con- 
secutive reactive ion etching in a parallel plate reactor. 
To find the most reliable parameters for the patterning 
process, they were optimized on 9 nm thick NbN samples 
structured into 130 nm wide stripes. As etching gases, 
pure SFg as well as SFg in combination with Ar or O2 
with varying flow ratios, gas pressures and etching power 
of the reactive ion etching process in conjunction with 
the electron beam exposure parameters have been inves- 
tigated. The criteria for the optimization were a mini- 
mal reduction of Tc in respect to the unstructured film 
and a minimal spread in the resistance around the nom- 
inal value. For processes fulfilling these criteria, they 
were further optimized for the maximal jc-values with 
the highest reproducibility. The optimal conditions were 
achieved for a SFgiAr mixture with 3:1 flow ratio, a total 
pressure of 90mTorr, an etching power of 100 W and 90s 
etching time. Those parameters were then used to struc- 
ture the 300 nm wide nanowire samples with the bends. 
In a consecutive step, the contact pads were made from 
the same film by means of photolithography and the same 
etching technique. The mean line thickness d = 10.4 nm 
was obtained by a step-profilometer as an average of sev- 
eral measurements at different locations along the line. 
The actual line-width in our samples varied between 250 
and 300 nm; it was measured for each sample with an 
accuracy better than 10 nm with a scanning electron mi- 
croscope. 



B. Results 

The superconducting transition temperatures of all 
studied samples were spread over the interval 13.1 - 13.6 
K. However we did not observe any correlation of the Tc 
values with either the bending angle or the radius of the 
connecting arcs. 

Figure [T] illustrates the effect of the bending angle on 
the measured critical current for samples from the batch 
#1. The critical current density jc{4:.2K) was calcu- 
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FIG. 1. (Color online) Averaged critical current densities at 
4.2 K for samples with different bending angles. Zero angle 
corresponds to the straight line. Inset shows the sample layout 
for the 60° angle. 



all samples have 90° bends, the critical current density 
increases with the increase in the inner radius r of the 
connecting arc. Kt r — 500 nm, the critical current of the 
line with bends amounted at almost 80% of the critical 
current in the straight line. 

We found that the bends affect the density of the crit- 
ical current differently at different temperatures. Fig. [3] 
shows temperature dependences of the critical current 
density for three samples from the batch #1 with bending 
angles 0°, 45° and 90°. These samples were selected for 
their smallest difference in the measured width to avoid 
additional impact of the sample width on jc (T)i^iil. The 
principal behaviour is, however, similar for all samples. 
In contrast to the significant difference in jc(4.2K) the 
jc values of all three samples practically coincide in the 
temperature range from the superconducting transition 
to approximately O.TTp. At lower temperatures, the 
jc values in lines with bends become lower than in the 
straight lines and the difference grows as the temperature 
decreases. 



III. DISCUSSION 
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FIG. 2. (Color online) Averaged critical current densities at 
4.2 K for samples with two 90° bends and different inner radii 
of the connecting arc. The inset shows the sample layout with 
the arc radius 500 nm. 



lated from the measured critical current at 4.2 K and the 
known width and thickness of each line. The straight 
lines without bends have jc(4.2K) = IS.BMA/cm^. For 
the samples with bends the critical current density grad- 
ually reduces as 9 increases. For samples with 90° bends, 
the critical current density is 8.1MA/cm^, which is ap- 
proximately 60% of the value for the straight line. 

Figure [2]shows the critical current densities j'(7(4.2K) 
for the samples from batch #2 with rounded 90-degree 
bends and different radii of connecting arcs. Although 



We compare the experimentally achieved critical cur- 
rent densities in the three samples from the batch #1 
with the depairing current density, which has been com- 
puted for each sample in the framework of the Ginzburg- 
Landau model as 



Jc(0) 



7^2 



21C(3)^/3 ep/Dlfci7b 



(1) 



with the temperature dependence proposed by Bardeeni. 
The current densities were further corrected for the ex- 
treme dirty limit according to the exact calculations of 
Kupriyanov and Lukichcv"^. We used the normal state 
resistivity p — 1.84]jilm of our lines at 20 K along with 
the electron diffusivity = 5 • 10~^m^s~^ and the en- 
ergy gap at zero temperature A(0) = 2kBTc; both are 
typical for 10 nm thick NbN film o^^i^^ . The results are 
shown in Fig. [5] where the ratio of the experimentally 
measured critical current density to the depairing critical 
current density jc / jc C^) is shown as function of temper- 
ature for the straight line (triangles) and for the line with 
two 90° bends (circles). It has been found in many ex- 
periments that even in straight lines the critical current 
reaches the depairing value only in the narrow interval 
below the superconducting transition. At lower tempera- 
tures in our straight sample experimental critical current 
amounts to «55% of the depairing value and this ratio 
remains almost temperature independent at T < 0.9 Tc". 
The line with the bends also carries an almost depairing 
supercurrent near the transition, but at low tempera- 
tures the ratio gradually decreases from more than 60% 
at T « 0.9 Tc to almost 40% at 0.3 Tp. 



4 




0.3 0.4 0.5 0.6 0.7 0.! 

T/T^ 



0.9 1.0 



FIG. 3. (Color online) Temperature dependences of the criti- 
cal current densities for three samples with different bending 
angles. Zero angle corresponds to the straight line. The er- 
ror of the critical current density is mainly defined by the 
accuracy of measurement of the width. The inset shows two 
typical current- voltage curves measured at 4.2 K for a straight 
sample (red) and a 90° sample (brown). 



We used the GL equations^"* to numerically calculate 
the critical current of the superconducting line with a 
rectangular cross-section and one sharp 90° bend. Simu- 
lations were made for the line width w ~ 60^gl(0) which 
is close to the experimental situation; in our films the 
Ginsburg-Landau coherence length at zero temperature 
^gl{0) is approximately 5nroi^. Since the widths of the 
studied lines in our experiment were all much smaller 
then Acfi'^", we neglected the screening effects. Although 
the Likharev's limit is fulfilled almost in the full tem- 
perature range except very close to Tc, the magnetic 
field generated by the applied transport current is not 
strong enough to overcome the edge barrier for vortex 
penetratio n -'^^1^^ . We therefore assume that the samples 
are initially in the Meissner state which becomes unsta- 
ble at some current which we define as the critical cur- 
rent. This is confirmed by an analysis of the critical 
current dependence on external magentic field applied 
perpendicular to the sample surface. It has been shown 
theoretically^^, that in thin and narrow {d ^ w Xcff) 
superconducting strips at small fields the critical current 
is 



Ic{H)=Ic{0)-—H 
•Acff 



(2) 



if the strip is in the Meissner state. Fig. [5] shows typical 
Ic{H)/ lc{0) dependences measured at 4.2 K for bridges 
with different widths. The bridges were made from a 
superconducting film with a thickness similar to the one 
from which the samples with bends were structured. The 
solid lines are linear fits of the experimental curves at 
i? — >■ 0. The slope of the curves increases with the width 
of the bridge. The inset in Fig. [5] shows the values of 
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FIG. 4. (Color online) Critical current in the line with two 
sharp 90° bends relative to the critical current in the straight 
line (filled diamonds) and the theoretical predictions of the 
GL (solid line) and London (dashed line) models. Open tri- 
angles and circles show the critical currents in the straight 
line and in the line with bends, respectively, relative to their 
depairing critical currents. 



the derivative at the linear part, which is proportional to 
(solid line) in good agreement with Eq. 2. The red 
point corresponds to the dIc{H)/dH value obtained on 
the straight bridge from batch #1. 

The reduction of the critical current in the framework 
of the London model was computed with Eq. 108 of Ref. 
13 describing the sharp bend. The solid and dashed lines 
in Fig. m depict the relative critical current (the ratio 
of critical currents in the line with bends to that in the 
straight line) obtained with the Ginzburg-Landau model 
and the London model (see Eq. 108 of Ref. 13), respec- 
tively. The theoretical calculations are limited to temper- 
atures T < 0.9 Tc since the London model is only valid 
for line-widths much larger than the coherence length 
(note that in the GL approach, the theoretical critical 
current approaches depairing current as T — >■ Tc). Al- 
though both models predict a decrease in the relative 
critical current at low temperatures, we found that the 
GL critical current is larger (by a factor of 1.3 for the 
given width) than it follows from the London model. We 
believe that the difference appears due to controversy 
in defining the energy barrier for single vortex entry in 
the framework of the London model. Indeed, there is a 
free parameter, which formally describes the minimal dis- 
tance at which vortex may approach the edge of a straight 
line. For consistency the minimal distance is defined by 
setting the barrier to zero for the current equal to the GL 
depairing current. Although slightly different in different 
publications, this distance appears comparable with the 
size of the vortex core. On the other hand, the barrier 
in the London model can be defined only assuming that 
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FIG. 5. (Color online) Dependence of normalized critical cur- 
rent Ic/IciQ) on an external magnetic field perpendicular to 
the sample surface for strips with different width as listed in 
the legend. The solid lines are linear fits of Ic{H) for H ^ 0. 
The inset shows the dependence of Ale /AH on width of su- 
perconducting strips. Symbols are experimental results, the 
solid line is a uP'-fit according to Eq. (2). The red point 
marks the value for the 300 nm wide sample from batch #1 
with 6* = 0°. 



this distance is larger than the radius of the core. Our 
calculations for narrower line [w — 20^gl(0)) show that 
the GL and London critical currents differ by the same 
factor «1.3 only at T < 0.6 T^; at larger temperatures 
this factor decreases and depends on the ratio w /^gl{T). 

Numerical and analytical results of the London 
modellS, which we discussed above, demonstrate that a 
bend becomes active in reducing the critical current when 
the line width is much larger than the coherence length. 
Therefore, decreasing the line width in SNSPD meander 
should improve the spectral range not only due to the 
decrease in the cross-section of the line but also due to 
weakening of the effect of bends on the critical current. 

Since even in the straight line the measured critical 
current is smaller than the depairing current, direct com- 
parison of jc{T)/j^{T) with the theoretical ratio is not 
wholly appropriate. To estimate the exclusive contri- 
bution of the bends, we use the ratio of the experi- 
mental critical current in the line with bends to that 
in the straight line (red diamonds in Fig. [3]). This ra- 
tio decreases with the decrease in the temperature and 
reaches approximately 0.65 at O.STc, whereas the theo- 
retical value are both below 0.5 (solid and dashed lines). 
Assuming that finite rounding may be present even in 
nominally sharp 90° bends and using the expressions 
for the reduction factor from Ref. 13, we have found 



negligibly small correction to the critical current for the 
rounding radius ^25 nm. This is the largest value of the 
bend radius which we have found when inspecting lines 
with sharp bends by the scanning electron microscope. 
Disagreement with the theory appears most likely due 
to the difference in the local homogeneity of two sam- 
ples (the origin of such non-homogeneities is not clear at 
the moment), which compete with the bends in reduc- 
ing the critical current. At this point we want to stress 
that jc{90°)/jc'(0°) values above one at temperatures 
T > 0.7 To stem from these competing effects between 
the local inhomogeneities and the rounding of the wires. 
The bends seem to dominate at T < 0.7 Tc where the 
relative critical current approaches theoretical values. 

To clarify the matter, further studies should be per- 
formed on more uniform (less granular) films from an- 
other material. For this case, it has been shown that 
for straight lines, the depairing critical current can be 
measured^Si. Further measurements of the effect of bends 
on such nanowires may allow for more consistent compar- 
ison between the theory and the experiment. 

IV. CONCLUSION 

In summary, we have shown that bends in supercon- 
ducting lines reduce the experimentally achievable crit- 
ical current and that the amount of the reduction im- 
posed by particular bend increases with the increase in 
the bending angle and the decrease in the rounding ra- 
dius of the bend. In spite of weaker reduction in the 
critical current predicted by the Ginsburg-Landau model 
as compared to the London model, we have experimen- 
tally found values which are less than expected from 
any considered model and have different temperature de- 
pendence. As the reason for this discrepancy we have 
proposed local non-homogeneities in the superconduct- 
ing lines. 

After submission of our paper, a very similar study 
of critical current in NbTiN films with 90° bends was 
published^^. The found reduction of the critical cur- 
rent was much smaller than the theoretical prediction 
of the London model^'^ and closer to the result given by 
the Ginzburg-Landau approach (if one assumes that the 
same correction factor wl.S is valid for wider films than 
ours). 
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